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ABSTRACT 


Twc-parameter  geowetric  systems  are  developed  for  noses 
and  tails  attached  to  parallel  middle  bodies  in  terms  of  inde¬ 
pendent  pol/Tiomials.  A  new  parameter  is  introduced--the  rate 
of  change  of  curvature  with  arc  length  at  the  juncture  with 
the  parallel  middle  body.  “Quadratic"  polynomials  are  con¬ 
sidered  for  bodies  of  revolution  and  "square  root"  polynomials 
for  symmetric  two-dimensional  bodies.  Permissible  ranges  of 
the  two  adjustable  parameters  are  obtained  for  selected 
geometrical  constraints  such  as  the  presence  of  inflection 
points. 


ADMINISTRATIVE  INFORMATION 


The  work  described  in  this  report  was  sponsored  by  the  Naval 
Ordnance  Systems  Command  (Code  035B)  and  was  funded  under  UR  123-01-03. 

INTRODUCTION 


The  geometrical  characteristics  of  streamlined  bodies  of  revolution 
and  symmetrical  two-dimensional  bodies  have  been  developed  in  a  previous 
report.^  There,  a  body  is  divided  into  a  forebody  and  an  afterbody  at  the 
maximum  section,  and  the  curvature  at  the  junction  is  one  of  two  adjustable 
parameters  controlling  the  shape.  These  forebodies  or  afterbodies  may  also 
be  attached  to  parallel  middle  bodies,  in  this  case,  it  is  considered  more 
desirable  in  hydrodynamic  applications  if  the  junction  curvature  is  made 
zero  to  match  that  of  the  parallel  middle  body  in  order  to  avoid  discon¬ 
tinuities  in  curvature.  The  use  of  a  junction  curvature  of  zeio,  however, 
eliminates  an  adjustable  parameter  and  leaves  only  one  adjustable  parameter 
for  the  polynomials  involved. 

It  is  the  purpose  of  this  report  to  rfcvurn  to  a  two-parameter  system 
for  forebodies  or  noses  and  afterbodies  or  tails  to  be  attached  to  parallel 
middle  bodies  in  order  to  have  a  more  extensive  group  of  shapos,  A  new 
adjustable  parameter  has  been  introduced,  namely,  the  rate  of  change  of 


^Granville,  P.  S.,  "Geometrical  Characteristics  of  Streamlined  Shapes," 
NSRDC  Report  2962  (^^r  1969);  Journal  of  Ship  Research,  Vol.  13,  No.  4 
(Dec  1969). 
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curvature  at  the  junction  with  the  parallel  middle  body.  The  requirement 
of  zero  curvature  here  becomes  an  additional  restraint. 

The  method  developed  in  Reference  1  is  followed  closely.  The  so- 
called  "quadratic"  and  "square  root"  polynomials  are  used  in  a  system  of 
independent  polynomials  controlled  by  the  two  adjustable  parameters.  The 
analytical  description  of  the  polynomials  is  obtained  by  a  factorial  method 
where  feasible.  Rounded,  pointed,  and  cusped  ends  are  considered.  Per¬ 
missible  ranges  of  the  adjustable  par.uneters  are  considered  in  terms  of 
geometric  criteria  such  as  the  presence  of  inflection  points.  Regions  of 
well-behaved  shapes  are  delimited  by  boundary  curves  obtained  by  an 
envelope  analysis.  The  least  squares  fit  of  graphically  or  analytically 
delineated  shapes  is  also  considered.  The  fullness  of  the  shapes  is  ob¬ 
tained  in  terms  of  a  prismatic  coefficient. 

GENERAL  ANALYSIS 

The  shapes  of  families  of  noses  or  tails  of  bodies  of  revolution 
and  of  two-dimensional  symmetrical  bodies  may  be  stated  functionally  as 

Y=frx;a,,  63  (1) 

■'  j  =  1,  2,  ... 

where  Y  is  the  radius  of  the  body  of  revolution  or  the  offset  of  the 
two  dimensional  body, 

X  is  the  axial  distance  of  the  body  measured  from  the  end  of  the 
nose  or  tail, 

are  the  parameters  to  be  varied  which  specify  the  family,  and 
Gj  are  the  boundary  conditions  or  restraints. 

The  analytical  analysis  is  more  useful  in  a  normalized  coordinate 
system  [x,  yj:  y  =  0  at  x  =  0,  and  y  =  1  at  x  =  1. 

If  D  is  the  diameter  or  maximum  thickness  of  the  parallel  middle 
body,  and  L  is  the  length  of  the  nose  or  tail,  then  the  normalized 
coordinates  become 


y  = 


2Y 

D 


(2) 


2 


and 


X 


(3) 


as  shown  in  Figure  1. 

To  achieve  "hydrodynamic  continuity"  as  contrasted  with  mathematical 
continuity,  it  is  considered  necessary  that  the  position,  slope,  and  curva¬ 
ture  match  at  the  junction  of  the  nose  or  tail  with  the  parallel  middle 
body.  In  this  case  the  slope  and  curvature  at  the  junction  are  zero. 

In  normalized  coordinates  the  contour  is  given  by 

i  =  1,  2,  ... 

y  =  f  [x;  a.,  B.]  (4) 

^  j  =  1.  2,  ... 

where  and  B^are  now  defined  in  normalized  coordinates.  In  this  study  two 
adjustable  parameters  and  02  are  to  be  considered  for  simplicity  of 
analysis. 

If  a  functional  form  like  that  of  a  polynomial  is  selected  as 

.n=n  ^ 

y  =  Z  a  x  =  P  [x]  (5) 

n=0 

a  resolution  into  linearly  independent  polynomials  may  be  obtained  like  that 
of  linearly  independent  vectors  multiplied  by  scalars,  such  as 

y  =  Z  f[a  ]  .  [x]  +  Z  f  [B.]  P„  .  [x]  (6) 

^  i  IlfX  .  J  II, J 


or  since  the  Bj  are  constant 


y  =  Z  f  [a.]  .  [x]  +  Q[x] 


whore 


Q[x]  =  Z  f  IBj]  Pjj^.  [x] 


(7) 


(8) 


3 


This  permits  the  effect  of  the  controllable  parameters  a.  to  be  obtained 
independently  of  each  other. 

The  independent  polynomials  may  be  determined  by  substituting,  one 
by  one,  conditions  and  6^  into  the  general  polynomial,  by  evaluating 
the  polynomial  coefficients  by  a  solution  of  the  resulting  simultaneous 
algebraic  equations,  and  then  by  a  gathering  of  terms  corresponding  to 
each  a^.  Another  method  is  to  use  the  factorial  properties  of  polynomials 
in  considering  conditions  and  6^,  which  is  illustrated  in  the  specific 
cases  to  follow. 

Not  all  combinations  of  give  desirable  shapes.  It  will  be 
of  value  to  analyze  the  possible  limitations  in  terms  of  simple  criteria: 

1.  Zero  condition- -y  =  0  for  0  x  £  1;  negative  values  of  y  would 
be  meaningless. 

2.  Unity  condition— y  =  1  for  0  £  x  1;  bulges  above  y  ■  1  are 
undesirable. 

3.  Maximum  or  minimum  condition- -dy/dx  =  0  for  0  £  x  £  1;  no  other 

maximum  or  minimum  is  to  be  permitted  than  at  x  =  1. 

2  2 

4.  Inflection-point  condition— d  y/dx  =  0  for  0  ^  x  £  1;  inflection 
points  are  undesirable  on  noses. 

For  example,  the  condition  for  zero  values  of  y 

y  [x;  Oj,  02]  =  0  0  X  1  (9) 

may  be  studied  as  follows. 

If  ttj  and  02  are  now  considered  as  variables,  and  x  is  considered 
as  an  adjustable  parameter,  a  line  may  be  defined  for  each  x-value.  An 
envelope  to  these  lines  may  be  developed  which  represents  the  boundary 
of  regions  for  values  of  and  values  of  y  =  0  for  different 

values  of  x.  The  envelope  condition  is  given  by 

■|x  ^  “r  “2^  ^  ° 
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From  Equations  (9)  and  (10) 


[X] 

(11) 

^2  M 

(12) 

A  plot  of  Oj  against  0,2  for  the  range  of  values  of  x,  0  £  x  1  gives  the 
envelope  curve.  Other  conditions  may  be  handled  in  a  similar  way. 

For  the  shapes  to  be  considered,  a  common  adjustable  parameter  to  be 
used  is  the  rate  of  change  of  curvature  with  arc  length  at  the  junction  of 
the  nose  or  tail  with  the  parallel  middle  body  at  x  =  1,  given  by 

El  -  (f)  ^ 

where  k  is  curvature,  and  £  is  arc  length. 

In  general 


k  = 


(14) 


and 


1  + 


(15) 


For 


X 


1, 


dy  _  d  y  _ 


dx 


dx" 


0 


and  then 


(16) 
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"QUADRATIC"  POLYNOMIAL  REPRESENTATION 


GENERAL 

The  functional  relation 

,  n=N 

y  =  I  a  X  (17) 

n=0 

is  called  the  "quadratic"  polynomial  for  uant  of  a  better  name.  It  is  very 
suitable  for  describing  bodies  of  revolution  for  which  volume  is  an  impor¬ 
tant  consideration,  since  it  represents  the  axial  distribution  of  the 
cross-sectional  area.  It  has  the  additional  advantage  of  providing  a  means 
of  accommodating  the  analytical  description  of  bodies  with  rounded  ends, 
something  the  ordinary  polynomial  cannot  do  because  of  the  requirement  of 
infinite  slope  at  x  =  0.  In  addition  to  bodies  with  rounded  ends,  the 
"quadratic"  polynomial  may  also  be  applied  to  bodies  with  pointed  ends  and 
cusped  ends. 

Although  any  number  of  adjustable  parameters  may  be  used,  the 
analysis  here  is  to  be  limited  to  two.  This  is  sufficiently  general  for 
describing  a  wide  variety  of  geometrical  shapes. 

ROUNDED  ENDS 

The  adjustable  parameters  to  be  considered  at  each  end  of  the 
curve  are 

Uj:  r  =  radius  of  curvature  at  x  =  0 


7 


02:  Kj  =  rate  of  change  of  curvature  at  x  =  1 


. . 


(19) 


The  boundary  conditions  6^  are 


S^:  X  =  0,  y  =  0 


X  =  1,  y  =  1 


^3-  ^  ^  ° 


(20) 


d^y 

Q:  x  =  l,-2-§  =  0 

^  dx^ 


Since  there  are  six  conditions  in  all,  n  =  5. 

The  and  6^  are  substituted  into  the  polynomial.  Differentiating 
Equation  (17)  successively  with  respect  to  y  gives 


2y  =  (aj  ^  2  a2X  +  3  ajX^  +  4  a^x^  +  5  a^x^)  ^ 


(21) 


and 


2 

2  =  (a^  +  2  a2X  +  3  a^x^  +  4  a^x^  +  5  a^x'*) 

dy 

+  (2  a2  +  6  a^x  +  12  a^x^  +  20  agX^)^-^^ 


(22) 


Since  a^  f  0,  dx/dy  =  0  at  x  =  0.  This  automatically  provides  a  rounded 
end.  Then 


a, :  a,  =  2r 

1  X 
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The  other  substitutions  yield 


6  Sj  24  +  60  =  2 


Sr  •o  =  '> 


h-  “O  *  “l  *  “2  *  “3  *  *4  *  “s  '  * 


a^  +  2  a2  +  3  aj  +  4  a^  ♦  5  ag  «  0 


8^:  2  a2  +  6  aj  +  12  a^  +  20  ag  =  0 

The  solution  of  Equations  (23)  by  determinants  shows  that  the  a's 
are  linear  functions  of  r  and'k^.  Hence  y  is  also  a  linear  function  of 
r  and  kj^  and  may  be  written  as 

y2  =  r  R[xj  ♦  kj  iCjtx]  +  Q[x]  (24) 

where  R[x],  Kj^[x],  and  Q[x]  are  also  polynomials  of  the  fifth  degree  in  x. 
It  is  possible  to  determine  R[x],  Kj[x],  and  Q[x]  by  first  solving  for  the 
a's  from  the  simultaneous  equations  and  then  regrouping  terms  applicable  to 
R[x],  K^[xl,  and  Q[xl.  Another  method  is  to  use  polynomials  as  follows. 

It  is  evident  that  the  relations  for  and  correspond  to 


“i- 


j3  ^  " 

“2=  fr y  ^  h 

dx 


Bj:  y^[0]  =  0 
^2'  1 
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(25) 


^3=  =  0 

64:  -^y^m  =  0 

Since  the  foregoing  apply  identically  to  r  and  because  of 
linearity,  it  is  further  evident  that 

(Xy  R^[0]  =  2,  Kj^[0]  =  Q^[0]  =  0 

a2:  Kj^^^tl]  «  2,  R^^^[l]  =  Q^^^[l]  =  0 
Bj:  R[0]  =  Kj[0]  =  Q[0]  =  0 
62:  Q[l]  *  1.  R[13  =  Kj[l]  =  0 

63:  R^[l]  =  Kj^Cl]  =  Q^[l]  =  0 

B4:  r"[1]  =  Kj“[l]  =  q”[1]  =  0 

where  R^  =  dR/dx,  R^^  =  d^dx^,  etc. 

Evaluation  of  R[x] 

Since  R[0]  =  R[l]  =  R^[l]  =  R^Ifi]  =  ^  ^ 

polynomial  of  the  fifth  degree,  R[x]  may  be  written  as 

R[x]  =  a  x  (x  -  I)**  (27) 

Since  R^[0]  =  2,  a  =  2. 

Then 

R[x]  =  2  X  (x  -  1)^  (28) 
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Evaluation  of  Kj[x] 


Since  K,[0]  *  Kj^[0]  »  Kj[l]  =  »  0,  k^[x]  may  be 

written  factorially  as 

Kj[x]  =  e  x^ 

Since  =  2,  B  =  1/3,  then 

KjW 

Evaluation  of  Q[x] 

Since  Q^[0]  =  Q^[l]  =  Q^^[l]  =  Q^^^[l]  =  0,  Q^[x]  may  be  written 
factorially  as 

Q*[x]  =  Y  X  (X  -  1)^  (31) 

Then  integrating  produces 

/  5  „  4  ,  2  \ 

Q[x]  =  Y  ^  ^  f^2) 

With  Q[0]  =  0  and  Q[l]  =  1,  C  =  0,  and  y  ■=  20.  Then 

Q[x3  =  -  (4  x^  -  15  +  20  X  -  10)  =  1  -  (x  -  I)"*  (4  x  +  1) 

(33) 

For  Rounded  Ends  in  Summary 

y^  =  r  Rtx]  +  kj  kj[x]  +  Q[x]  (24) 

with 

R[x]  =  2  X  (x  -  1)"^  (28) 


(X  -  !)• 


(X  -  !)• 


(29) 


(30) 
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(30) 


I 


Kilx]  Cx  -  1)^ 

Q[xJ  =  1  -  (X  -  I)'*  (4  X  ♦  1)  (33) 

PERMISSIBLE  RANGES  OF  PARA^^ETERS  r  AND 

Zero  Condition 

-  f[x;  r,  kj]  *  0  0  <  x  <  1  (34) 

The  envelope  in  r  and  k^  with  x  as  the  variable  parameter  is  given  by 

f‘-||=0  (3S) 

The  two  envelope  conditions.  Equations  (34)  and  (35),  provide  two 
simultaneous  equations  in  r  and  kj  which  are  solved  by  the  Cramer  rule  to 
give  r(x]  and  kj[x]: 


x^  (3  x^  -  10  X  t-  10) 
2  (X  -  l)'^ 


k  -i-iii 


5  x^  + 


(X  -  1) 


10  X  - 
3 


101 


(36) 


(3?) 


The  envelope  curve  is  shown  in  Figure  2.  Since  the  tangent  lines  represent¬ 
ing  values  of  r  and  k^  for  y  =  0  are  outside  the  envelope  curve,  the  inside 
region  contains  values  of  r  and  kj  which  do  not  have  values  of  y  =  0  (except 
at  X  =  0). 

Unity  Condition 

The  unity  condition  is  that 

y^  =  f  [x;  r,  kj]  =  1  0  1.  x  S  ^ 


i 


4 

i 

I 

$ 

I 

] 

i 

1 


L. 


12 
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The  envelope  in  r  and  with  x  as  the  variable  parameter  is  given  by 


|j  Cf  -  I)  =  0  (39) 

The  two  envelope  conditions.  Equations  (38)  and  (39),  provide  two 
simultaneous  equations  in  r  and  kj,  which  are  solved  by  the  Cramer  rule  to 
give 


3  X  +  2 

2  X 


(40) 


^1  = 


3  (X  -  1)‘ 


(41) 


The  envelope  curve  is  shown  in  Figure  2.  Desirable  values  of  r  and  k^, 
that  is,  without  bulges,  are  on  the  "inside  curved"  side  of  the  envelope 
curve. 


Mar.icum  or  Minimum  Condition 

The  maximum  or  minimum  condition  is  given  by 

=  f^  =  0  (42) 

The  envelope  curve  in  r  and  k^  with  x  as  the  variable  parameter  is  given 
by 


=  0  (43) 

The  envelope  curve  is  shown  in  Figure  2.  A  better  understanding  of 
the  envelope  curve  is  developed  in  Figure  3.  Each  point  on  the  envelope 
curve  represents  a  tangent,  giving  the  locus  of  values  of  r  and  kj^  which 
provide  a  maximum  or  minimum  at  each  value  of  x  other  than  the  maximum  at 
X  =  1  which  prevails  at  all  times.  Two  such  loci  are  represented.  Their 
point  of  intersection  provides  a  value  of  r  and  k^,  representing  maxima  or 
mir.ima  at  two  values  of  x.  Evidently  from  any  point  in  the  region  outside 
the  envelope  curve,  two  tangents  may  be  drawn  to  the  envelope  curve.  Thus 
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0  1C  20  30  40  60 


figure  2  -  "Quadratic"  Polynomial:  Rounded  End,  Pennissible  Range 

o£  Parameters  r  and 


Figure  3  -  "Quadratic"  Polynomial:  Rounded  End,  Delineation  of  Regions 
by  Envelope  Curve  for  Maximum  or  Minimum  Condition 
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the  region  outside  the  envelope  curve  represents  values  of  r  and  giving 
two  maxima  or  minima.  The  region  inside  the  envelope  cur\'e  provides  no 
maximum  or  minimum.  Finally  there  is  only  one  maximum  or  minimum  specified 
by  the  envelope  curve. 

The  two  envelope  conditions,  Equations  (42)  and  (43),  provide  two 
simultaneous  equations  in  r  and  k^,  which  are  solved  by  the  Cramer  rule 
to  give  r[x]  and  kj[xJ  as 


r 


15  X 


10  X  -  5  X  +  1 


(44) 


W  -  50  (X  -  1)^ 

i^l  -  2 

^  10  X  -  5  X  +  1 


(45) 


Inflection- Point  Condition 

The  inflection-point  condition  is  given  by 


For  y^  =  f[x] 


II  I 

2  f  r  -  r  =  0 


(46) 


and  the  envelope  condition 


f^^^  =  0  (47) 

The  two  conditions  provide  two  simultaneous  equations  in  r  and  kj^  in  terms 
of  X,  Since  the  boundary  condition  leads  to  a  quadratic  relation,  the  Cramer 
rule  does  not  apply.  For  specified  values  of  x  the  two  simultaneous  equa¬ 
tions  may  be  numerically  solved  by  direct  substitution  of  one  equation 
into  the  other.  The  results  are  shown  in  Figures  2  and  4.  The  results  for 
the  rounded  end  are  summarized  in  Table  1. 
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Inflection-Point  Condition 


TABLE  1  -  SUKWARY  OF  "QUADRATIC"  POLYNOMIALS  FOR  ROUNDED  END 

Polynoiiltl  Equttlon: 

•  r  R[x)  ♦  k,  K,tx3  ♦  Qlx]  0  <  x  <  1 

R[x3  •  2  X  (x  -  D* 

K,Cx]  -^x^  (X  -  1)3 


Q[x]  -  1  -  (x  -  D*  «  X  ♦  1) 
En«t1ope  Equations: 

Zero  Condition: 

.2  ,2 


.  x^  (3  X*  -  10  X  ♦  10)  ,  r  .  3  (x3  -  5  x^  ♦  10  X  ♦ 

2  U  -  D*  '  ’  (•. 


JOI 


Unity  Condition: 

r.l2LL  ;  i,  .  3 

2  X 

HexiMi  or  Klnlotfs  Condition: 

.2 


r  • 


IS  x‘  .  j  .  30  (x  -  D' 

- - 

10  x^  -  5  X  ♦  1  10  x^  -  5  X  ♦  1 


Inflection-Point  Condition:  IlMrlctl  c«1cj1*t1on 


POINTED  ENDS 


The  adjustable  parameters  are 


a, :  s  =  slope  at  x  =  0 


a2‘-  kj  =  rate  of  change  of  curvature  at  x  =  1 


(48) 


(49) 


(20) 


The  boundary  conditions  6^  are  the  same  as  those  given  in  Equation 


6^:  X  =  0,  y  =  0 
6^:  X  =  1,  y  =  1 


(50) 


Since  the  "quadratic"  polynomial  as  previously  used  giver  infinite  slope  at 
X  =  0,  an  additional  condition,  indicated  as  follows,  is  necessary  to  give 
controlled  slopes  at  x  =  0.  Hence,  the  degree  of  the  polynomial  becomes 
six. 

Forcj 


dy  a,  +  2  a-  X  +  . . .  +  n  a  x^  ^ 
_ _ 1 _ 2 _ n _ 

dx  “  2  y 


(51) 
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Since  y  =  0  at  X  =  0,  dy/dx  «,  unless  =  0.  For  a^  =  0,  dy/dx  is 
indeterminate  at  x  =  0.  Then  by  the  L'Hopital  rule 


X 


L 

00 


dy 


(52) 


or 


s^  = 


Oj  then  requires  that 


a-  =  0 
1. 


^2  =  ^ 


The  other  substitutions  yield 


(53) 


6  a3 

+  24  a^  +  60  ag  + 

^0  = 

0 

^2  = 

^3  = 

^  * 

2  a2  +  3  Sj  +  4  a 

^4^ 

2  a.2 

+  6  a_  +  12  a.  +  i 
3  4 

(54) 


2  2 
y  is  then  a  linear  function  of  s  and  or 


y^  =  s^  S[x.l  +  Kj[x]  +  Q[x] 


(55) 


The  analysis  proceeds  as  before,  and  the  results  are  shown  in 
Table  2  and  Figures  5  and  6. 
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TABLE  2  -  SlJ^WARY  OF  "QUADRATIC"  POLYNOMIALS  FOR  POINTED  END 
Polynomial  Equation: 

S[x]  +  K^Cx]  +  QCx]  0  <  x  <  1 

where 

S[x]  «  x^  (x  -  1)^ 

K^Cx]  »  jx^  (x  -  1)^ 

Q[x]  «  1  -  (x  -  1)^  (10  x^  +  4  X  +  1) 

Envelope  Equations: 

Zero  Condition: 

2  ,  3  x^  (2  x^  -  6  X  4  S)  .  r  ^  6  (2  x^  »  9  x^^  15  x  -  M 
'  (x-1)«  (x-1)5 

Unity  Condition: 

2  ,  3  (2  x^  2  X  1)  .  I  .  6  (x  -  1)^  (2  X  +  1) 
s  *  2  *  *^1  ~  3 

X  X 

Maximum  or  Minimum  Condition: 

.2  _  30  x^  .  t  .  60  (x  -  1)2 

5x^-4x+l  *  5x^-4x  +  l 

Inflection-Point  Condition:  Numerical  calculation 
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40 


Figure  5  -  "Quadratic"  Polynomial: 

Permissible  Range  of  Parameters  s^  and  kj 


Pointed  tnd, 
2 


Figure  6  -  "Quadratic"  Polynomial:  Pointed  End, 
Inflection-Point  Condition 
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In  analyzing  the  inflection-point  condition,  a  problem  arises  for 

A 

X  =  0,  where  an  indeterminate  condition  exists.  By  the  L'Hopital  rule 


-  24  s  -  2  +  120 


X  =  0 


The  boundary  curve  at  x  =  0  is  then 


-  12  s  -  kj  +  60  =  0 


CUSPhD  ENDS 


The  adjustable  parameters  are 


a,:  k  =  curvature  at  x  =  0 
1  o 


X  =  0 


a^:  kj  =  rate  of  change  of  curvature  at  x  =  1 


X  =  1 


The  boundary  conditions  Bj  are  the  same  as  those  of  Equation  (20), 
except  for  the  additional  condition  of  zero  slope  at  x  =  0. 

6j:  x  =  0,  y  =  0 


®2'  X  =  0-  ^  “ 


Bji  X  =  1,  y  =  1 
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®4*  ^  "ax " 


Be*  X  =  1,  — ^  =0  (6 

As  will  be  shown  the  cusped  end  requires  two  additional  conditions  which 
makes  n  =  8. 


dy  a,  +  2  a-  X  +  ...  +  n  a  X 


n  -  1 


"2*  dx 


2  Caj  X  .  ...  ♦  X  ) 


.  0  at  X  .  0,  aj  .  0 


dy  2  a-  +  3  a,  X  +  ...  +  n  a  "  ^ 

^  _ 2 _ 3 _ n _ 

dx  "  ^2 

2  (a2  *  aj  X  +  ...  +  a^^  X  '  ) 


(*).. 


(s) 


=  0  ,  a,  =  0 


X  =  0 


For 


d^y  6  X  +  ...  +  n(n  -  1)  ~  ^  -  2 


.2 


dx 


2  y 


(65) 


Then 


d^y  2[6a_+. . .n(n-l)a  x”  ^](a_+...+a  x”  (  3a_+...+na  ) 

_ j _ n  -*5 _ n _ \  5 _ n _ / 


dx 


.  I  1/3  n-8/3 

4  (a  x  •'^...+a  X 
3  n 


3/2 


(66) 


2 

Equation  (66)  gives  d^’y/dx  -»•  “>  at  x  =  0. 
To  prevent  this,  let  a^  =  0. 

Now 


d^y  2[12a^+. .  .+n(n-l)a^x""^]  (a^+. .  .+aj^x""'^)- ^  4a^+. .  .+naj^x""^  ) 


dx" 


At  X  =  0 


4(a^+...+aj^x"‘'*  ) 


3/2 


(67) 


=  2  a 

,2  4 

dx 


(68) 


4  a.  =  k  " 
4  o 


(69) 


The  other  substitutions  yield 


O2:  6  a^  +  . . .  +  n  (n  -  1)  (n  -  2)  a^^  =  2  kj^ 

B,:  a.  +  . . .  +  a  =1 
3  4  n 
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B*:  4a.  +...+na  =0 

4  4  n 

65:  12  a^  +  ...  +  n  (n  -  1)  =  0  (70) 

2  2  ■“ 
y  is  then  a  linear  function  of  k  and  k, 

^  o  1 

y^  =  K^[x]  +  kj  Kj[x]  +  Q[x]  (71) 

The  analysis  proceeds  as  before,  and  the  results  are  shovm  in  Table  3 
and  Figures  7  and  8. 

"SQUARE  ROOT"  POLYNOMIAL  REPRESENTATION 

GENERAL 


The  functional  relation 

y  =  a,/2  £  a.  -  (72) 

is  to  be  called  the  "square  root"  polynoniial  for  want  of  a  better  name. 

It  is  suitable  for  describing  two-dimensional  shapes  with  rounded  ends, 
since  the  square-root  term  gives  infinite  slope  at  x  =  0.  Of  course, 
without  the  square-root  term  an  ordinary  polynomial  remains. 

The  same  analysis  procedure  used  for  the  quadratic  polynomial  is  to 
be  applied  where  possible  to  the  square-root  polynomial  for  the  same  cases: 
rounded  ends,  pointed  ends,  and  cusped  ends. 

ROUNDED  ENDS 

The  adjustable  parameters  are  the  same  as  those  for  the  quadratic 
polynomial 


a^:  r  =  radius  of  curvature  at  x  =  0 
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TABLE  3  -  SUMMARY  OF  "QUADRATIC"  POLYNOMIALS  FOR  CUSPED  END 


Pciyncmial  Equation: 

•  ^0^  Kj,[x]  +  K^[x]  +  Q[x]  0  <  X  <  1 

where 

Kq[x]  =  }  x^  (x  -  1)^ 

K^[x]  =  y (x  -  1)^ 

Q[x3  *  X®  (35  x^  -  120  x^  +  140  x  -  56) 

Envelope  Equations: 

Zero  Condition: 

2  _  4  x^  (15  x^  -  40  X  +  28)  .  _  12  (5  x^  -  20  x^  +  28  x  -  14) 

Unity  Condition: 

.  2  ^  4  (15  x^  -f  20  x^  18  x^  12  X  5) 

0  4 

X 

j;  .  12  (x  -  1)^  (5  x^  5  x^  3  X  1) 

1  ^ 

Maximum  or  Minimum  Condition: 

k  2  =  840  x^  .  r  ,  840  (x  -  1)^ 

°  14  x^  -  16  X  +  5  ’  ^  14  x^  -  16  x  +  5 


Inflection-Point  Condition:  Numerical  calculation 


Figure  8  -  "Quadratic"  Polynomial:  Cusped  End 
Inflection- Point  Condition 


^2‘  ~  change  of  curvature  at  x  =  1 


:  ^  ] 

[dxy 


(74) 


X  =  1 


The  boundary  conditions  6^  are  the  same  as  those  for  the  quadratic 


polynomial 


gj:  X  =  0,  y  =  0 


$2'  X  =  1,  y  =  1 


e  :  X  =  0 

dx 


(75) 


Since  there  are  six  conditions  in  all,  n  =  4. 

The  and  g^  are  substituted  into  the  polynomial.  For  a^: 
Differentiating  Equation  (72)  with  respect  to  y  gives 


or 


1  = 


4 


"l/2  ^ 


-1/2  . 


n  -  K  dx 

n  a  X  )  -r- 

n  ^  dy 


dy 


1/2 


la  +  a  + 

2  ^1/2 


. . .  +  n  a  X 
n 


2n  -  1/2 


(76) 


(77) 


At  X  =  0,  dx/dy  =  0  which  ensures  a  rounded  end.  Differentiating  Equation 
(76)  with  respect  to  y  gives 


+  n  a  x~ 
_ n 


2n 


+  n  a  X 

_ a 


2n  -  1/2 


4  ®l/2 


-  n  (n  -  1)  X 


2n  -  1/2 


(78) 
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At  X  =  0,  r  =  1/2  a' 


1/2 


h/2  = 


(79) 


or 


2’ 

8  ^1/2 

k  • 

a  =  0 

1 

0 

‘2  = 

"1/2^ 

1 

’3* 

2  ^1/2 

-  1  a 

’4' 

4  1, 

-  *^1 


‘1  ^  “2  "  “3  "  **4 


=  0 


(80) 


The  presence  of  the  square-root  term  prevents  the  use  of  the 
factorial  analysis. 

The  solution  as  simultaneous  equations  in  the  a's  produces 


y  =  /Tr  R[x3  +  kj  K^[x]  +  Q[x] 


(81) 


with 


R[x]  =  (5  -  21  x^  +  35  X  -  35) 


(82) 


KlM  =-|-  (x  -  1)^ 


(83) 


Q[x]  =  1  -  (X  -  1)"^  (84). 

The  permissible  ranges  of  parameters  /2  r  and  are  studied  as 
before,  and  the  results  are  shown  in  Table  4  and  Figures  9  and  10. 
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TABLE  4  -  SUMMARY  OF  "SQUARE  ROOT"  POLYNOMIALS  FOR  ROUNDED  END 


Polynomial  Equation: 

y  ■  Jl~T  R[x]  +  [x]  ♦  Q[x]  0  1  x  i  1 


where 


R[x]  -  /T  +  *g  (5  x'^  -  21  x'^  +  35  X  -  35) 


K^Cx]  ••I'X  (x  -  1)- 


Q[x]  «  1  -  (x  -  1) 
Envelope  Equations: 

Zero  Condition: 


28  X  -  4  -  (3  x'^  -  14  X  +  35) 


r  .  24  (  7  x^  -  20  x^  *  18  X  -  4)  -  3  x^^^  (x^  -JO  x~  ■*■  59  x^  -  112  x  *  70) 


(x  -  1)^  [28  X  -  4  .  (3  x^  -  14  X  +  35)] 


Unity  Condition: 


_8jx-jj^x^ 


^  x^  [28  X  -  4  -  x^'*^  (3  x^  -  14  X  +  35)] 


n 


3  [x  -  1)  [56  x^  ■*■  8  X  -  x^^  (x^  -  7  x^  35  X  •*•  35)] 


x^  [28  X  -  4  -  x^^^  (3  x^  -  14  X  +  35)] 


Maximum  or  Minimum  Condition: 


T 


(x-l)Vx3/g 


28  x*-  -  17  X 


-  48  (x 

^  /IQ  w3  ~  £9 


-  62  ♦  91  X  -  35) 


j;  ,_6  J[x  -  1L[28  x^4_-jl^/LL3_xI-  14  x,-^  35)1 


’  28  x^  -  17  X  +  1  ~  x^/^  (18  x^  -  62  x^  +  91  x  -  35) 

Inflection- Point  Condition; 


^~r  ’ 


96  (x  -  1)^  x^^^ 


14  x^  -  15  X  ♦  3  +  2  x^^^  (18  -  40  X  ♦  21) 


'  14  x*^  -  15  X  ♦  3  +  2  x’'^  (18  x^  -  40  X  + 


21) 
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POINTED  ENDS 


The  ordinary  polynomial  is  utilized,  namely 


y  = 


n=n 

Z 

n=0 


a 


n 


X 


n 


(85) 


The  adjustable  parameters  ct^  are  the  same  as  those  for  the  quadratic 
polynomial 


s  =  slope  at  x  =  0 


■(«).. 


(86) 


X  =  0 


a^:  kj  =  rate  of  change  of  curvature  at  x  =  1 


^1  = 


d  y 


(87) 


dx 


X  =  1 


The  boundary  conditions  3^  are  the  same  as  those  for  the  quadratic 


polynomial 


3j:  X  =  0, 
62=  X  =  1, 

V  X  = 

3^:  X  =  1, 

Since  there  are  six  conditions  in  all. 


y  =  0 
y  =  1 


^  = 
dx 


d"y 

dx^ 


n  =  5. 


(88) 
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Substitution  of  and  0^  into  the  polynomial  produces 


II 

s 

“2  = 

6  a3 

+  24  a.  +  60  ac  = 
4  5 

01 : 

0 

^2  = 

^1  "  ^2  "  ^3  "  ^4 

^1  " 

2  a2  +  3  a^  +  4  a 

^4  = 

2^2 

+  6  a,  +  12  a.  + 

3  4 

(89) 


V  is  a  linear  function  of  s  and  k,,  or 

1 


y  =  s  S[x]  +  kj  Kj[x]  +  Q[x] 


(90) 


The  analysis  proceeds  as  before,  and  the  results  are  shown  in  Table  5 
and  Figures  11  and  12. 

CUSPED  ENDS 

Th?  ordinary  polynomial  is  utilized 


n=n 

y  =  E  a  x‘ 
n  n 
n=0 


n 


(91) 


The  adjustable  parameters  are  the  same  as  those  for  the  quadratic 
polynomial 


a, :  k  =  curvature  at  x  =  0 
1  0 


■  (21 . . 


(92) 
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TABLE  5  -  SUMMARY  OF  ORDINARY  POLYNOMIALS  FOR  POINTED  END 
Polynomial  Equation: 

y  =  s  S[x]  +  K^[x]  +  Q[x]  0  £  x  £  1 

where 

S[x]  =  X  (x  -  1)"^ 

(X  -  1)3 

Q[x]  =  1  -  (x  -  1)^  (4  X  +  1) 

Envelope  Equations: 


Zero  Condition: 


Figure  11  -  Ordinary  Polynomial:  Pointed  End, 
Permissible  Range  of  Payaneters  s  and 


Figure  12  -  Ordinary  Polynomial:  Poisited  End,  Inflect ion- Point  Condition 


kj  =  rate  of  change  of  curvature  at  x  =  1 


(93) 


X  =  1 


The  boundary  conditions  6^  are  the  same  as  those  for  the  quadratic 


polynomial 


Bji  X  =  0,  y  =  0 


Qy  X  =  1,  y  =  1 


(94) 


^4*  *  "  ° 


d^y 

Br*  X  =  1,  ~  0 

^  dx^ 


Since  there  are  seven  conditions  in  all,  n  =  6. 

Substitution  of  and  Bj  into  the  polynomial  produces 


“r 


a~:  6  a-  +  24  a.  +  60  ar-  +  120  a,  =  k, 

2  3  4  5  6  1 


6r  "o^"" 


62=  =  0 


(95) 


Bj:  ag  +  a^  +  a^  +  +  a^  +  a^  +  a^  =  1 


B;,:  a  +  2  a-  +  3  a,  +  4  a.  +  5  a^  +  6  a,  =  0 
4  1  2  3  4  5  6 


B^:  2  a2  +  6  a^  +  12  +  20  a^  +  30  a^  =  0 
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y  is  a  linear  function  of  and  k,  in 
'  o  1 


y  =  kjj  K^[x]  +  kj  Kj[x]  +  Q[x] 


(96) 


The  analysis  proceeds  as  before  and  the  results  are  shown  in  Table  6  and 
Figures  13  and  14. 

LEAST- SQUARES  FIT 

Given  shapes  either  in  graphical  form  or  in  analytical  form  may  be 
fitted  to  the  polynomials  of  this  report  by  a  least-squares  fit.  If 


z 


for  quadratic  polynomials,  or 


z  =  y 

for  square-root  polynomials;  then,  in  general 

z  =  d  D[x]  +  kj  Kj[x]  +  Q[x]  (97) 

where  d  and  D[x]  refer  to  the  appropriate  type  of  shape.  For  example, 
d  =  r,  and  D[x]  =  R[x]  for  the  rounded  nose. 

1  2 

In  general,  a  least-squares  fit  requires  that  (r.  -  z^)  dx 
be  minimized  where  2j^[x]  represents  the  body  shape  to  be  fitted.  Conse¬ 
quently,  by  differentiating  with  respect  to  the  coefficients  to  be  deter¬ 
mined,  there  result  two  simultaneous  equations 

d  /J  dx  +  kj  /J  D  iCj  dx  +  /J  D  Q  dx  -  /J  Zj  D  dx  =  0  (98) 

d  /J  D  Kj  dx  +  kj  dx  +  /j  Kj  Q  dx  -  /J  z^  dx  =  0  (99) 

The  values  of  d  and  are  then  determined  by  the  Cramer  rule.  Values  of 
the  indicated  integrals  for  the  various  polynomials  may  be  found  in  Table  7. 
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TABLE  6  -  SU^WARY  OF  ORDINARY  POLYNOMIALS  FOR  CUSPED  END 
Polynomial  Equation: 

y  =•  ko  Kq[x]  +  [x]  +  Q[x]  0  <  X  <  1 

where 

Kj,[x]  -  J  x2  (X  -  1)^ 

KiCx]  =  6  (X  -  1)^ 

QCx]  =  1  -  (x  -  1)^  (10  x^  +  4  X  +  1) 

Envelope  Equations: 

Zero  Condition: 

.  6  (2  •  6  »  *  5)  .  J  ,  12  (2  x^  •  9  x^  *  15  x 

(x-1)«  •  ’  (X  -  1)3 

Unity  Condition: 

.  _  6  (2  x^  2  X  1)  .  Z  .  12  Oc  -  1)^  (2x^-1) 

o'  x2  ’  1  ■  x3 

Maximum  or  Minimum  Condition: 

k  ,  ,  W  xf  .  .  .  J  .  JIO  ijL2Li)L 

0  5x^-4x  +  l  ’  5x'^-4x  +  l 

Inflection-Point  Condition: 

1(  .  60  x^  (10  x^  -  10  X  3) 

°  50  x^  -  80  x^  +  45  x^  -  10  X  +  1 

k  -  120  (X  -  (10  x^  -  5x^1) 

^  50  X  -  80  x^  +  45  x^  -  10  X  +  1 
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Figure  14  -  Ordinary  Polynomial:  Cusped  End,  Inflection-Point  Condition 
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SQUARE  ROOT"  "QUADRATIC 


TABLE  7  -  INTEGRATED  VALUES  FOR  LEAST- SQUARES  FIT 


Polynomial 


Rounded 


B 

B 


Pointed  S 


X 

CM 

a 

O 

Jo  “ 

f  D  Q  dx 

Jo 

Jo 

Jo 

4 

1 

14 

1 

43 

495 

■  1980 

495 

20,790 

1 

1 

202 

1 

1 

529 

6435 

■  3U;55ff 

mM 

■  366,360 
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PRISMATIC  COEFFICIENT 


The  fullness  of  a  shape  is  given  by  the  prismatic  coefficient,  which 
is  the  ratio  of  the  volume  of  the  body  to  the  volume  of  a  prism  having  the 
maximum  cross-sectional  area  and  the  length  of  the  body. 

The  prismatic  coefficient  Cp  for  the  shapes  of  this  report  is  then 

Cp  =  /J  y^  dx  (100) 


for  bodies  of  revolution,  and 


C 


P 


dx 


(101) 


for  two-dimensional  bodies. 

For  the  polynomials  of  this  report 

Cp  =  d  /J  D  dx  +  kj  \  dx  /J  Q  dx  (102) 


Consequently  for  bodies  of  revolution  described  by  quadratic  polynomials 
rounded  end 


C 


P 


(103) 


pointed  end 


cusped  end 


2520 


1512 


(104) 


(105) 
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ard  for  two-dimensional  bodies  described  by  square-root  polynomials 


rounded  end 


pointed  end 


cusped  end 


Lines  of  constant  C 


5 


C 


P 


2 

3 


C 


P 


k 

0 

210 


+ 


7 


are  plotted  in  Figures  4,  6,  8,  10,  12,  and  14. 


(106) 


(107) 


(108) 
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